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We applied the mean-field approach to a columnar bundle assembled by the parallel arrangement 
of stiff polyelectrolyte rods in a salt bath. The electrostatic potential can be divided into two 
regions: inside the bundle for condensed counter-ions, and outside the bundle for free small ions. To 
determine the distribution of condensed counter-ions inside the bundle, we use a local self-consistent 
condition that depends on the charge density, the electrostatic potential, and the net polarization. 
The results showed that, upon bundle formation, the electric charge of polyelectrolytes, even those 
inside the bundle, tend to survive in an inhomogeneous manner, and thus their width remains finite 
under thermal equilibrium because of the long-range effect of charge instability. 



I. INTRODUCTION 

For electrically charged soft materials including bio- 
logical polyelectrolytes such as DNA and actin, it is im- 
portant to understand their behavior as a many-body 
Coulomb system. A well-knovifn model under a mean 
field approximation is the Poisson-Boltzmann equation 
(PB eq.), which has often been used to describe the dis- 
tribution of small ions around a stiff polyelectrolyte rod 
in a good solvent A typical example is negatively 
charged DNA with a persistence length of ~ 50 nm. The 
PB eq. allows one to estimate an ion distribution not 
only for a short rigid fragment of DNA, but also for a 
segment of a long coiled DNA shorter than the persis- 
tence length. However, we have to be very careful when 
interpreting the results of the PB for polyelectrolytes in 
a poor solvent. 

DNA condensation is induced by various condensing 
agents such as multivalent cations, neutral poly- 

mers (polyethylene glycol, etc), and cationic surfactants. 
The experimental observation of the aggregation of short 
oligomeric DNA suggests that their local structure is 
an aligned ordered packing f^. Furthermore, this ap- 
plies not only to condensed DNA, which consists of the 
multiple-assembly of short fragments, but also to com- 
pact DNA, which consists of a single long chain. Upon 
the addition of condensing agents, long DNA discretely 
folds from a coil to a compact state at the level of a 
single molecule. This DNA compaction can be charac- 
terized as a first-order phase transition under Landau's 
argument @ . The typical morphology of compact DNA 
is a toroid, the local structure of which can be regarded as 
an ordered bundle based on the parallel arrangement of 
its segments, when we note the cross-section of a toroid 
cut with the plane including the rotational symmetric 
axis. 
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The transition to a parallel packed arrangement should 
lead to a drastic change in the distribution of small ions. 
The condensation of counter ions inside the packed state, 
which accompanies both the bundling of short DNAs and 
the folding of a long DNA chain, should be enhanced to 
diminish coulombic instability. In fact, the experimental 
results support this significant change Q • The condensed 
counter-ions in compact DNA have a notable contribu- 
tion, but this is not appropriately estimated by the PB 
eq. This implies that a new approach that is different 
than the PB eq. is needed for this case. 

In this article, we adopt a different approach from 
the PB eq. for a bundle that consists of the parallel- 
arrangement packing of stiff polyelectrolyte rods in the 
presence of salt. The morphology also includes a col- 
umn with a finite diameter, since the experimentally ob- 
served size is finite. For example, the bundle assembled 
by stiff rod-like polyelectrolytes such as F-actins has a 
finite width Q. Also, in the case of a long single DNA 
chain, it has been found that the coil and compact states 
coexist in a single molecule under some conditions [9]. 
In addition, at high concentrations of long DNA, a stiff 
thick bundle with a persistence length on the order of /im 
by multiple-assembly has been observed At present, 
it is difficult to experimentally determine whether these 
phenomena regarding a finite size arise from the free en- 
ergy minimum in the equilibrium condition or a kinetic 
effect. In this study, we obtain the minimum point at a 
finite diameter on the free energy profile, and discuss the 
mechanism. 



II. COLUMNAR BUNDLE ASSEMBLED BY 
THE PARALLEL ARRANGEMENT OF STIFF 
POLYELECTROLYTE RODS IN A SALT BATH 

Imagine a columnar bundle that is formed by the par- 
allel arrangement of stiff polyelectrolyte rods, as shown in 
Fig 1. A stiff polyelectrolyte rod carries charge —Q{< 0) 
and has a line density of charge —q/d [q is the unit 
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FIG. 2: (Color Online) Fraction of condensed counter-ions per 
polyelectrolyte 9{r) as a function of the radial distance from 
the central axis of the bundle r. This profile is obtained from 
the differential equation (Eq. ([6])) at various values of electro- 
neutrality at the central axis 6{r — 0). (a) 9a., (b) 9h,(£) 9^, 
(d) 6»d, (e) fo. (fa < < < < 9,). See reference UM for 
these values. The parameters are d = 0.17 nm, Ib ~ 0.71 nm, 
a = 1 nm, v = 1 nm^, A — 0.1 nm. 



FIG. 1: (Color Online) Schematic representation of a colum- 
nar bundle assembled by the parallel arrangement of stiff 
polyelectrolyte rods. Monovalent counter-ions (monovalent 
cations) condense into the bundle, while counter-ions and co- 
ions are in the bulk. The remaining electric charge of the 
bundle is screened by the surrounding small ions in bulk, r 
is the distance from the central axis of the columnar bundle, 
and R is the bundle radius. 



charge), which means that the length of the rods is 
L (= Qd). A bundle composed of these rods is dissolved 
in the presence of a monovalent salt (1+,!—), the con- 
centration of which is c. The intrinsic dielectric con- 
stant in solvent (water) is eo- Also, we introduce the 
Bjerrum length Ib = / A-Ke^k-QT = q^/S/iTreo (fee is the 
Boltzmann constant, T is the thermal temperature, and 
/3= IAbT). 

In this system, monovalent counter-ions can not induce 
the bundling of polyelectrolytes independently without 
any additional attraction. As an additional attractive 
force, we consider non-electrostatic forces such as the de- 
pletion force by a neutral polymer, etc. This situation 
corresponds to DNA compaction induced by polyethy- 
lene glycol (PEG). 

To formulate this system, we use the concept of 
the Oosawa-Manning condensation theory (OM the- 
ory), which has been applied to a dispersed polyelec- 
trolyte 12] by dividing the potential into two regions: 
inside and outside the polyelectrolyte. This theory is su- 
perior to the PB eq. not only because it is easy to under- 
stand, but also because it corresponds to the experimen- 
tal results. Analogously, we divide the system into two 
regions, with small ions (counter-ions and co-ions) out- 
side the bundle, and condensed counter-ions inside the 
bundle. 

Here, we consider the self-consistent local condition for 



the interior ion distribution, by relating D{r), E{r) and 
P{r) in the following relation: 

D{r) = eoE{r) + P{r), (1) 

where D{r) is the electric displacement field, E{r) is the 
electric field, P{r) is the polarization, and r is the radial 
coordinate perpendicular at the central axis of the bun- 
dle, as shown in Figure [TJ Hereafter, we regard D{r), 
E{r) and P{r) as the "mean" and "continuous" field on 
a coarse-grain scale. The cross-sectional density of poly- 
electrolyte rods in the bundle are also uniform, and the 
cross-sectional area per polyelectrolyte is v/d. The con- 
densed small ions into the bundle are monovalent cations, 
the fraction of which is defined as 9{r) per polyelec- 
trolyte rod. In the following discussion, we consider that 
< e{r) < 1. 



1. Electric field: E{r) 

In the equilibrium state, the chemical potential is 
equivalent everywhere inside the bundle, as follows: 

A^i^JcW+KmW^ constant. (2) 

where is the chemical potential of the mean elec- 

tric field inside the bundle and fi]^g^{r) represents the 
contribution of translational entropy. Here, we neglect 
the higher-order multipole composition. The electro- 
static chemical potential is written as fJ-"l^{r) — q(f>{r). 
The chemical potential for the translational entropy of 
a monovalent counter-ion is also given by /iira(^) = 
log {9 {r)/cv). Therefore, the derivative of the to- 
tal chemical potential with respect to r is d{^]^^{r)) + 
^i^a(r))/5r = 0. Thus, this leads to 
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qmr) = a.(log^) = dr{\og0{r)), (3) 
cv 

where the mean electric field E{r) = —d<j){r)/dr. 



2. Polarization: P(r) 

The degree of polarization should be closely related 
to the excluded volume of polyelectrolytes and counter- 
ions. Counter-ions should condense on the surface of each 
polyelectrolyte rod in the bundle. We may assume that 
the separation of the two charges of the dipole is the 
closest-approach distance. In this case, the polarization 
density is considered to be 



P{r)^-qXeir)/v, (A > 0), 



(4) 



where A is the effective mean approach between the poly- 
electrolyte and counter-ions. 
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3. Electric displacement field: D{r) 



From the definition of 0(r), the charge density is given 



divl?(r) = -q{l - e{r))/v. 



4- Local condition inside the bundle 



(5) 



div I?(r) = eodiv£'(r) + divP(r) is given by the diver- 
gence in Eq (H]). If we include Eqs. ©, (g]) and (O, we 
obtain the local condition of the ion distribution in the 
bundle as follows: 
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lBr{l-e{r)) = —dr{rdr{logeir))} + lBdri-rX0ir)). 
47r 

(6) 



5. Boundary condition at the central axis (r = 0) 

In the above equation, there are two conditions at the 
central axis (r = 0) and the bundle surface {r — R). The 
boundary condition at the central axis can be obtained 
under the assumption that the displacement field inside 
is continuous everywhere. In this case, the displacement 
field is zero at the central axis. Therefore, D{r = 0) = 
leads to 



D{r = 0) = eoE{0) + P{0) 

kBT dr9{r)\r^Q qXeiO) 



9(0) 



0. (7) 



Diameter of bundle (nm) 

FIG. 3: The plot of the free energy jSfk/kL as a function of the 
bundle diameter 2R (= 2y^n7/7rd) with changes in the value 
of the surface tension 7 at A = 0.1 nm. fk/k corresponds to 
the free energy per stiff polyelectrolyte rod in the bundle with 
width R. The parameters used in the calculation are 7 = 0, 
0.120, 0.197, 0.200, and 0.203 nm"^ at a salt concentration 
c = 0.3 M. The other parameters are the same as those in 
Fig.H 



where Dr{r) is the radial displacement field [ij]. Note 
that D{r) is represented by E{r) and P{r) (Eqs. ([3]), 
©)• 

In addition, the boundary at r = i? will be naturally 
associated with the relation to the total surviving charge 
of the bundle, which is screened by the linearized PB-type 
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ion atmosphere. Note that the total surviving charge in 
the bundle will be determined by the free energy mini- 
mum as discussed later. 

Figure [5] shows examples of solutions at various pro- 
portions of condensed counter-ion at the central axis, 
9{r = 0). The proportion of condensed counter- ion inside 
the bundle greatly decreases far from the central axis. 



A. Free energy of the bundle with a finite width 

In this section, we consider the bundle with a finite 
width, using the above local condition. First, let us con- 
sider the coexistence of gas and liquid in a finite- volume 
system that contains a finite number of particle, to em- 
phasize the essential points regarding the finite-width 
phenomenon of a bundle. Gas and liquid coexist in the 
presence of some attraction such as van der Waals force. 
In this situation, the finite cluster sizes arise from the 
penalty of translational entropy when the number of gas 
molecules in the system decreases. However, the mecha- 
nism in the observed bundle with a finite width is clearly 
different. 

An entropic penalty is incurred along with bundle for- 
mation by the condensation of small ions into the bundle, 
and by the decrease in the degrees of freedom on poly- 
electrolyte rods. However, a finite distribution of the 
bundle width is observed under the experimental condi- 
tions, where sufficient small ions are dissolved to act as 
a particle bath [8] , and where there is a high concentra- 
tion of rods 15]. Therefore, to essentially understand the 
finite- width distribution of the bundle as an electrostatic 
effect, we must consider the mechanism, which should 
not be attributed to the effect of the finite number of 
polyelectrolyte rods or small ions. 

In this system, the total free energy Fe-tot is written 

as 



B-cnt 



Fn- 



B-int- 



(8) 



-FB-cnt corresponds to the degrees of freedom of the poly- 
electrolyte ( translational and rotational entropy for stiff 
rods, or the elastic free energy in the case of a semi- 
flexible chain). Fe-int is the free energy of the effective 
interaction part between rods. As mentioned above, we 
only examine the contribution that is essential for the 
finite- width effect Fe-int- This can be divided into dif- 
ferent contributions, as follows: 

^B-int = ^ nkfk 
k 

= ^ f^fe (/fe,clc-in + /fe,scr + //c.tra + /fe,atr) , (9) 



/3/fc.cie-i„ = L / 2iTrdT-——Dr{r) {Dr{r) - P{r)) , 



/3/fe,tra = L / 27rrdr {e{r)/v) log {e{r)/ecv), 
Jo 

/3/fe,atr = 27ri?L7, 



(10) 
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where fk is the free energy of the bundle formed by k 
stiff polyelectrolyte rods, and is the number of bun- 
dles labeled by k (total number of rods is krik), and 
kv — ■nR^d. fk,eie-in IS the electrostatic free energy in- 
side the polyelectrolyte bundle and /fc,scr is the contri- 
bution of the ion atmosphere around the bundle. The 
inverse Debye length is k = ^/Stt/bc and K^{x) are the 
modified Bessel functions of order v. /fc^tra is the contri- 
bution from the translational entropy of the monovalent 
countcr-ion. fk, a.tr is the additional attraction between 
polyelectrolyte rods, which is not attributable to small 
ions (e.g., depletion force by neutral polymer). The ion 
atmosphere around the bundle is calculated as the lin- 
earized PB-type, like that around a single highly charged 
polyelectrolyte [16]. Also, 7 is the surface tension (di- 
vided by /crT). We do not need to consider the volume 
part by attractive interaction, since that part of /fc^atr/^ 
is independent of k under the assumption of pair-wise 
attraction. When we apply the ion distribution obtained 
from Eq. (|6]), we minimize the free energy fk/k, which 
is the total free energy of one bundle fk divided by the 
constitutive number of polyelectrolyte rods. 



III. RESULTS AND DISCUSSION 

Figure [3] shows the profile for the free energy of one 
bundle divided by the constitutive number of polyelec- 
trolyte rods and by its length (3fk/kL, as a function 
of i?, at various values of 7. The free energy fk/k at 
each width R (or k) is minimized by changing the degree 
of electro-neutrality at the central axis 6{r — 0). Note 
that the free energy is divided by k-QT{^ P^^) in this 
plot. As shown in Fig. [3l when the surface tension 7 
is greater, minimum points with a broad concave region 
appear. This indicates that bundles with a finite width 
are distributed. Moreover, these shift to a larger radius 
R with an increase in the surface tension 7. 

Also, Fig. |4](A) shows the free energy plot [3fk/kL as 
a function of the diameter 2R at various concentrations 
of salt c. Figure [11(B) shows the diagram of the bundle 
width at a minimum free energy /3fk/kL as a function 
of the salt concentration c. These results indicate that 
the minimum point shifts to a larger value when the salt 
concentration increases. 

Note that, if the translational and rotational entropy 
of single rods, or the elastic free energy of a semi-flexible 
chain arc added to the free energy fk/k, double min- 
imum points should appear in the free energy profile. 
This means that finite-width bundles and single rods co- 
exist in bulk, or folded and unfolded states coexist in a 
single chain. In fact, in the case of long DNA, single- 
chain segregation has been experimentally observed in 
the presence of a neutral polymer such as PEG ^ . 
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FIG. 4: Dependence of the bundle diameter on the salt con- 
centration. (A) The plot of free energy Pfk/kL as a function 
of the bundle diameter 2R with changes in the salt concentra- 
tion c at a fixed surface tension 7 = 0.203 nm~^. The values 
for the salt concentration are c = 0.260, 0.280, 0.290, 0.295, 
and 0.300 M. (B) Diagram of bundle diameters at minimum 
free energies as a function of the salt concentration. The other 
parameters used in the calculation are the same as those in 
Figs. H H 



Single-chain segregation has also been confirmed in the 
presence of multivalent cations Moreover, finite 



widths of a multiple aggregated assembly have been ob- 
served |15|]. These widths are close to those in our results, 
although they must be distinguished from PEG with re- 
gard to electrostatic effects. This should not be a sur- 
prise, even with the different condensing agents, since one 
of the first factors is the excluded volume effects, which 
are comparable in both cases. Furthermore, in the pres- 
ence of a multivalent cation, the growth trend of width 
upon DNA compaction is consistent with those in their 
multiple assembly. Similarly, this feature is expected to 
be valid for a PEG system, and our model satisfies this 
point. 

How does the bundle have a finite width? From the so- 
lution of the differential equation ^ , a charge in the inte- 
rior of the bundle is not completely neutralized by mono- 
valent counter-ions. This remaining charge is almost de- 
termined by the competition between electrostatic free 
energy and translational entropy [9]. In fact, the free 
energy minimum of the bundle shifts to a greater width 
with an increase in the salt concentration as shown in 
Fig.m This implies that the condensation of monovalent 
counter-ions into the bundle incurs a penalty in transla- 
tional entropy, and this entropic penalty decreases at a 
higher salt concentration. 

However, this contribution alone is not sufficient to 
discuss the finite-width effect. The electrostatic penalty, 
including the translational entropy of small ions, must in- 
crease slowly. If the electrostatic penalty is immediately 
moderated with bundle growth before the decrease in the 
surface penalty, the bundle will continue to grow toward 
infinite width. To clarify this point, we consider a sim- 
plified situation under the assumption that the bundle is 
uniformly neutralized by the ratio 9 without polarization. 
Electro-neutralization is almost completely determined 
by competition between the following two dominant con- 
tributions: 



M^T^l = L{7rR'9/v)logi0/ecv), 



(11) 
(12) 



where f^^^ is the electrostatic free energy in the bundle, 
and /^^a is contribution in the translational entropy 
of the counter-ion. Note that we use the electric field 
E{r) — —rq{l — 9)/2eQV obtained from Gauss's law in 
the calculation of f^^°- Here we expand it with respect 
to (5 = 1 — 0, since a thick bundle is almost neutralized. 
This minimization gives S ~ log {1 / cv) / {ttIbR'^ /2v + 1) ~ 
2v\og{l/cv)/{TTlBR'^) at i? > ^/v/Ib- By substitu- 



tion, we have Pif^^" + fk^J/kL 



log {ecv)/d 



3{\og{cv))^v/TTdlBR^- The second term represents the 
free energy penalty with bundle growth. On the other 
hand, the penalty of the surface energy is /3/fc,atr/fc ^ 
j/R. These sums + /^tTa + //c,atr)/fc can have a 

maximum point, but not a minimum point in the free en- 
ergy profile. This indicates that the bundle grows toward 
infinite width. 
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Therefore, for the finite distribution of bundle width, a 
further effect of long-range instability should be required. 
Our local condition in Eq. (jG]) leads to non-uniform neu- 
tralization by the condensed counter-ions in the bundle, 
which gives rise to slow-growth instability. Moreover, the 
polarization P{r) might reduce the electrostatic free en- 
ergy loss in the integrand of Eq. (jlOp more than that loss 
in the case of invariant dielectric constant (P(r) = 0). 

Here, we emphasize the significance of a many-body 
effect on bundle formation. As a model, we suppose a 
column morphology even in a thin bundle. Under this 
assumption, the surface penalty is underestimated, or at- 
traction is overestimated because the number of neigh- 
boring rods is overestimated. This implies that, even 
if the two-body interaction is repulsive, many-body in- 
teraction could give an gain in free energy to form the 
bundle. Although the finite-width bundle formed by stiff 
rods in the presence of multivalent cations has been ar- 
gued to be due to a kinetic effect in the case of attractive 
two-body interaction between the parallel rods [llj, at 
the same time bundle formation should be examined be- 
fore the two-body interaction becomes attractive, since 
many-body interaction has advantageous effects by de- 
creasing the surviving charge to reduce electrostatic re- 
pulsion, and by increasing the number of neighboring 
rods with bundle growth. 

Before we conclude, we address the applicable scope 
of the present model. The electric field is determined 
only by the translational entropy of condensed counter- 
ions, as seen in Eq. ([3]). Therefore, this model can be 
used if the correlation by condensed counter-ions is not 
significant. One possible example is chromatin fiber, if 
the quadrupolar interaction of the nucleosome structure 
is attractive and decoupled from condensed counter-ions. 
Although our proposed approach should include a cer- 
tain level of contribution that is neglected in the simple 



PB eq., this is similarly invalid, e.g., in the presence of 
a multivalent cation, where ion-exchange between mono- 
valent and multivalent cations is expected to accompany 
the transition. However, if the expression of the uncon- 
sidered local correlation is known and is added as a new 
term to the local condition in Eq. the coarse-grain 
electric field is given by the number of condensed counter- 
ions. Thus, we hope that our approach can be expanded 
to various cases. 



IV. CONCLUSION 

In conclusion, we propose the application of the mean 
field approach for a bundle with a parallel arrangement 
of stiff polyelectrolyte rods. There is an inhomogeneously 
surviving charge of polyelectrolyte rods in the bundle in- 
terior, which is given by the local conditions which de- 
pend on the charge density, the electrostatic potential 
and the net polarization. The bundle width exhibits a 
finite distribution due to the long-range effect of electro- 
static instability. Furthermore, the stable finite width 
shifts to a larger value with an increase in the salt con- 
centration. 
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